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Abstract 

The construction of a consistent theory of quantum gravity is a problem in theoretical physics that 
has so far defied all attempts at resolution. One ansatz to try to obtain a non-trivial quantum theory 
proceeds via a discretization of space-time and the Einstein action. I review here three major areas 
of research: gauge-theoretic approaches, both in a path-integral and a Hamiltonian formulation, 
quantum Regge calculus, and the method of dynamical triangulations, confining attention to work 
that is strictly four-dimensional, strictly discrete, and strictly quantum in nature. 



1 



1 Introduction 



With the absence of a satisfactory quantum theory of gravity, a major puzzle in theoretical physics 
remains still unsolved. This article contains an overview and a comprehensive bibliography of past 
efforts to define a consistent theory of quantum gravity in four dimensions via an intermediate 
discretization. I will discuss only models with some concrete implementation of the dynamics of 
Einstein's theory, Lagrangian or Hamiltonian. 

One way of tackling the quantization problem non-perturbatively is to use discrete methods, 
in analogy with quantum field theories on a flat background. They of course have to be modified 
in the case of gravity, where the metric of space-time becomes itself a dynamical variable. The 
diffcomorphism invariance of the classical continuum theory is no obstacle in principle to the 
introduction of a discretization, as has been demonstrated by the success of such methods in 
describing two-dimensional Euclidean gravity in a path-integral approach. 

The approaches I will describe below have yielded a variety of interesting results on the discrete 
simulation methods themselves, as well as on the geometric properties of typical configurations, 
the role of gauge invariance, the phase structure, the inclusion of matter, to name just a few. They 
have not so far been successful in providing convincing evidence for the existence of a non-trivial 
four-dimensional quantum theory of gravity (neither, of course, have other methods). There are 
considerable technical difficulties in performing analytical computations and setting up sufficiently 
big numerical simulations, ft is a theoretical challenge to come up with more realistic models that 
may lead to a truly interacting gravitational theory. 

I have divided the discrete approaches into three categories, which 1 will discuss in the order of 
their chronological appearance: the gauge-theoretic formulations using connection variables, based 
on first-order descriptions of Einstein gravity, and two metric formulations, the quantum Regge 
calculus program and the more recent method of dynamical triangulations, which use simplicial 
instead of hypercubic lattice discretizations. They mainly deal with pure gravity, with the possible 
inclusion of a cosmological constant term and higher-order derivative terms, and in some instances, 
matter-coupling. 

This review covers about 200 papers, the earliest of which appeared in 1979 (not counting the 
1964 article by Leutwyler |l45| , who used a lattice approximation to perform a gravitational sum 
over histories a la Feynman). Much of the motivation for the discrete investigations comes of course 
from other sources, which for the most part I do not cite explicitly. The reader is referred to the 
individual articles for a more complete bibliography of the background material. 



2 Gauge- theoretic discretizations of gravity 
2.1 Lagrangian treatment. Introduction 

This area of research was inspired by the success of non-perturbative lattice methods in treating 
non-abelian gauge theories |173| . To apply some of their techniques, gravity has to be brought into 
a gauge-theoretic, first-order form, with the pure-gravity Lagrangian 

S[A,e] = / eAeAR[A], (f) 

J M 

where the so(3, f )-valued spin-connection Ajf (with curvature R) and the vierbein are con- 
sidered as independent variables. The important feature of (|l|) is that Aj t J is a gauge potential, 
and that the action - in addition to its diffeomorphism-invariance - is invariant under local frame 
rotations. Variation with respect to leads to the metricity condition, 

^\A} = d [A]+ A ij e i} = ^ ( 2 ) 

which can be solved to yield the unique torsion- free spin connection A = A[e] compatible with 
the e^j. There are some obvious differences with usual gauge theories: i) the action (|I|) is linear 
instead of quadratic in the curvature two-form Rj£, of A, and ii) it contains additional fields e^. 
Substituting the solution to (||) into the action, one obtains S[e] — J M d 4 x (det e) R, where R 
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denotes the four-dimensional curvature scalar. This expression coincides with the usual Einstein 
action J d 4 x y/detg R only for det e > 0. 

Most of the lattice gauge formulations I will discuss below share some common features. The 
lattice geometry is hypercubic, defining a natural global coordinate system for labelling the lattice 
sites and edges. The gauge group G is 50(3, 1) or its "Euclideanized" form 50(4), or a larger group 
containing it as a subgroup or via a contraction limit. Local curvature terms are represented by 
(the traces of) G- valued Wilson holonomies Un around lattice plaquettes. The vierbeins are either 
considered as additional fields or identified with part of the connection variables. The symmetry 
group of the lattice Lagrangian is a subgroup of the gauge group G, and does not contain any 
translation generators that appear when G is the Poincare group. 

When discretizing conformal gravity (where G = 50(5, 1)) or higher-derivative gravity in first- 
order form, the metricity condition on the connection has to be imposed by hand. This leads to 
technical complications in the evaluation of the functional integral. 

The diffeomorphism invariance of the continuum theory is broken on the lattice; only the local 
gauge invariances can be preserved exactly. The reparametrization invariance re-emerges only at 
the linearized level, i.e. when considering small perturbations about flat space. 



2.2 Smolin's lattice model 



The first gauge-theoretic model for lattice gravity is due to Smolin [186|, based on the continuum 



formulation of MacDowell and Mansouri [157 1 , with de Sitter gauge group 0(3,2) or 0(4, 1), and 
a Lagrangian of i? 2 -type, 

'xe^euKLRliR^, (3) 
where the components of R^ v are related to those of the usual curvature tensor by 

R" = R"±^A-4 e l)- ( 4 ) 

Although the underlying gauge potentials are o(3, 2)- or o(4, l)-valued, the action is only in- 
variant under the 6-dimcnsional subgroup of Lorentz transformations. The theory contains a 
dimcnsionful parameter /. The gauge potentials associated with the internal 5-direction are iden- 
tified with the frame fields e^, and the action can be decomposed into the usual Einstein-term (Q) 
plus a cosmological constant term with A ~ -k and a topological R A i?-term. Smolin analyzed its 
lattice discretization, and found both a weak- and a strong-coupling phase, with respect to the di- 
mensionless coupling constant g ~ . He performed a weak-coupling expansion about flat space 
and rederived the usual propagator. In the strong-coupling regime he found massive excitations 
and a confining property for spinors. 



2.3 Numerical implementation 



Caracciolo and Pelissetto ]63|, |6J, |65[ |66[ |67J performed a numerical investigation of the phase 
structure of Smolin's model. Using the compact group 50(5) and its associated Haar measure, their 
findings confirmed the two-phase structure: a strong-coupling phase with a confining property and 
presence of exponential clustering, and a weak-coupling phase dominated by a class of topological 
configurations, with vanishing vierbein. However, their Monte Carlo data (on 4 4 and 8 4 -lattices 
with periodic boundary conditions) indicated strongly that the transition was first-order, even if 
the measure was generalized by a factor of | det e\ p , p G [0, 150] E7|. 



2.4 Other gauge formulations 



Das et al pS lattice-discretized an 5p(4)-invariant Lagrangian due to West [194 , leading to a 



functional form S[U] = yTr (([/„ - ±Tr U n ) 2 ), where U n denotes a sum of double-plaquette 
holonomies based at the vertex n. Details of the functional integration were not spelled out. The 
square-root form of the Lagrangian does not make it very amenable to numerical investigations 
(see also the comments in |139|). 
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Mannion and Taylor [159 suggested a discretization of the (J e A e A i?)-Lagrangian without 
cosmological term, with the vierbeins e treated as extra variables, and the compactified Lorcntz 
group SO (4:) 

Kaku |l38| has proposed a lattice version of conformal gravity (see also [137, 13E[), based on 
the group 0(4,2), the main hope being that unitarity could be demonstrated non-perturbatively. 
As usual, a metricity constraint on the connection has to be imposed by hand. 

A lattice formulation of higher-derivative gravity, containing fourth-order i? 2 -terms was given 
by Tomboulis [18S], The continuum theory is renormalizable and asymptotically free, but has 
problems with unitarity. The motivation for this work was again the hope of realizing unitarity 
in a lattice setting. The square-root form of the Lagrangian is similar to that of Das et al. It 
is 0(4)-invariant and supposedly satisfies reflection positivity. Again, the form of the Lagrangian 
and the measure (containing a 5-function of the no-torsion constraint) is rather complicated and 
has not been used for a further non-perturbative analysis. 

Kondo [ 142 1 employed the same framework as Mannion and Taylor, but introduced an explicit 
symmetrization of the Lagrangian. He claimed that the cluster expansion goes through just as in 
lattice Yang-Mills theory, leading to a positive mass gap. 



2.5 Proving reflection positivity 

Various properties of gauge gravity models were analyzed by Menotti and Pelissetto in a series of 
papers in the mid-eighties. They first studied a discrete 0(4)-version of the Lagrangian (|j) Jl65| , 
including the 0(4)-Haar measure and a general local 0(4)-invariant measure for the vierbein fields, 
and showed that reflection-positivity holds only for a restricted class of functions. Furthermore, 
expanding about flat space, and after appropriate gauge-fixing, they discovered a doubling phe- 
nomenon similar to that found for chiral fermions in lattice gauge theory, a behaviour that also 
persists for different gauge-fixings. One finds the same mode doubling also for a flat-background 



expansion of conformal gravity [ 167 1 . In the same paper, they gave a unified treatment of Poincare, 
de Sitter and conformal gravity, and showed that reflection positivity for 0(4)-gravity (as well as 
for the two other gauge groups) holds exactly and for general functions only provided a signature 
factor sign(det e) is included in the Lagrangian. 

To ensure the convergence of the functional integration, one has to introduce a damping factor 



for the vierbeins in the measure, both for Poincare and conformal gravity [167]. An extension of 



the results of [167| to supergravity with the super-Poincare group is also possible |62)j. One can 



prove reflection positivity and finds a matching gravitino doubling in the perturbative expansion. 



2.6 The measure 

One can derive non-trivial conditions on the lattice measure by imposing the Slavnov- Taylor iden- 
tities in a perturbative lattice calculation. Starting from the general form of the measure for 
Poincare gravity, 

[] det(e p (n),e,(n),e p (n),e (T (n)) Ar / 16 n/(^W) J] de^n) J] dU L (n, n + »), (5) 

n^ivpcr nfM n,fi>0 n,/.i>0 

Menotti and Pelissetto |l66|[ performed a one-loop calculation and found that the Slavnov- Taylor 
identity can be satisfied for a particular choice of N, /i, /a (momentum expansion coefficients of /), 
and of the cosmological constant A. The solution still depends on a real parameter £ (related to a 
residual non-invariance under rotations) . This makes it difficult to draw any immediate conclusions 
on the structure of the full, non-perturbative measure. 



2.7 Assorted topics 

Caselle et al |70| proposed a lattice action that is genuinely Poincare-invariant, at the price of intro- 
ducing additional lattice "coordinate variables" . They also suggested a compact 0(5)-formulation 
which reduces to the Poincare form in the limit as the length of some preferred 0(5)-vector is 
taken to infinity, as well as a super-version involving the graded Poincare group. The same authors 
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in p9| put forward an argument for why doubling should appear in general gravity plus matter 
systems. 

Reisenberger [175 has recently suggested a gauge-theoretic path integral based on the Plcbanski 
action for Euclidean gravity. He discretizes the theory on a simplicial or hypercubic lattice with 
group- and algebra- valued fields. A metricity constraint needs to be imposed on the basic spin-1 
fields, which it turns out is difficult to treat exactly. 



2.8 Summary 

The gauge-theoretic Lagrangian lattice approaches are afflicted by a number of technical difficulties. 
Reflection positivity can be shown for some of the models, but generally requires the inclusion of 
a factor sign(det e) in the Lagrangian. Obtaining qualitative non-perturbative information about 
the phase structure requires a non-trivial measure input. The complicated functional form of the 
Lagrangian and the metricity condition that has to be imposed via a Lagrange multiplier and 
the corresponding functional integration do not make conformal and higher-derivative theories 
attractive candidates for numerical simulations. The compact version of Smolin's de Sitter gravity 
is still the simplest model, but its numerical investigation did not yield interesting results. 



2.9 Hamiltonian treatment. Introduction 

Relatively little work has been done on discretized Hamiltonian fomulations of gravity. This 
can in part be understood from the fact that the numerical methods available for lattice gauge 
field theories rely mostly on the Euclidean path-integral description. Unfortunately the relation 
between the Lagrangian and Hamiltonian quantizations for generally covariant theories without a 
fixed background is far from clear. 

In the usual metric formulation, the complicated non-polynomial form of the Hamiltonian 
constraint has been a long-standing problem. In this framework, neither the functional form of 
the quantum representations nor the nature of the quantization problems suggest that a discrete 
approach might yield any advantages. This situation has improved with the introduction of new 
Hamiltonian gauge-theoretic variables by Ashtekar |^2|, ^3) . 

By a Hamiltonian lattice approach one usually means a formulation in which the time variable 
is left continuous, and only the spatial 3-slices are discretized. In continuum gravity, the 3+1 
decomposition leads to the (non-Lie) Dirac algebra of the three-dimensional diffeomorphism gener- 
ators and the Hamiltonian constraint, associated with the deformation of three-surfaces imbedded 
in four-space. One usually requires this algebra to be realized in the quantum theory, without 
anomalous terms, for a set of self-adjoint quantum constraint operators, for some factor-ordering. 

Since a discretization of space-time breaks the diffeomorphism invariance, there is no reason to 
expect the Dirac algebra to be preserved in any discrete approach, even classically. This raises the 
question of whether and in what form part of the diffeomorphism symmetry can still be realized 
at the discrete level. Using gauge-theoretic variables, one can maintain the exact local gauge 
invariance with respect to the internal degrees of freedom, but there is no analogous procedure for 
treating the coordinate invariance. 



2.10 Hamiltonian lattice gravity 

There is a gauge-theoretic Hamiltonian version of gravity defined on a cubic lattice, which in many 
aspects resembles the Lagrangian gauge formulations described earlier. It also is virtually the only 
discrete Hamiltonian formulation in which some progress has been achieved in the quantization 



(see also [155 for a recent review) 



Renteln and Smolin [ 180 1 were the first to set up a continuous-time lattice discretization along 
the lines of Hamiltonian lattice gauge theory. Their basic configuration variables are the link 
holonomies U(l) of the spatial Ashtekar connection along the edges. The lattice analogues of 
the canonically conjugate pairs (A l a (x) , Ef (x)) are the link variables (U (I) a B , Piipj) > with Poisson 
brackets 

{U A B (n,a),U c D (m,b)} = 0, (6) 
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{Pi(n, a), U A C (m 7 b)} = -\- <W#aS nA B U B C (n, a), (7) 
{p l (n 1 a),p j (m, b)} = 5 nm S & ^eij k Pk(n,a), (8) 

with the S'f7(2)-generators satisfying [Ti,Tj] = 2ejj/-Tfe. Lattice links I = {n,d) are labelled by 
a vertex n and a lattice direction a. These relations go over to the usual continuum brackets 
in the limit as the lattice spacing a is taken to zero. In this scheme, they wrote down discrete 
analogues of the seven polynomial first-class constraints, and also attempted to interpret the action 
of the discretized diffeomorphism and Hamiltonian constraints in terms of their geometric action 
on lattice Wilson loop states. 



2.11 The measure 

There is a natural measure for the quantum theory, given by the product over all lattice edges 
of the Haar measures dg. However, since the Ashtekar connections A are complex- valued, the 
gauge group is the non-compact group 5*0(3, (D) = SL(2, C), and the gauge-invariant Wilson loop 
functions are not square-integrable. For the alternative formulation in terms of real SU(2)-variables 
(see below), these problems are not present. An alternative heat kernel measure dv for holomorphic 
5i(2, C) holonomies on the lattice was used in |14£, |99|. 



2.12 The constraint algebra 

This line of research was continued by Renteln [1TS[] , who proved that for a particular factor- 
ordering (all momenta to the left), the subalgebra of the discretized diffeomorphism constraints 
(smeared by lapse functions TV), 

Y,V utt {N,n) := ]>> latt (n,a)Tr(lV^M), (9) 

n n 

closes in the limit as the lattice spacing is taken to zero. This calculation was later extended to a 
variety of different symmetrizations for the lattice ope rato r and to an arbitrary factor-ordering of 
the form aTr (Ut)p+ (1 — a)pTr (Ut), with < a < 1 |156|] . Again, one does not find any quantum 
anomalies. It would be highly desirable to extend this result to commutators involving also the 
discretized Hamiltonian constraint and to find the explicit functional form of the anomalies, if 
there were any. 



2.13 Solutions to the Wheeler-DeWitt equation 

Part of solving the canonical quantum theo ry is to determine the states annihilated by the Hamil- 
tonian constraint H. It was shown by Loll |l46|] that solutions exist in the Renteln-Smolin formu- 
lation, where H{n) = b e yfc Tr (U(n, U^r^p^n, a)pj(n, b). They are given by multiple, non- 
intersecting Polyakov loops (the lattice is assumed to have compact topology T 3 ). Such solutions 
are trivial in the sense that they correspond to quantum states "without volume" . The difficulties 
one encounters when trying to find other solutions is illustrated by the explicit calculations for the 
1 x 1 x 1-lattice in l!4£ 



The search for solutions was continued by Ezawa []99j (see also 100 for an extensive review), 
who used a symmetrized form of the Hamiltonian. His solutions depend on multiple plaquette 
loops (Ua) k , where a single lattice plaquette □ is traversed by the loop k times. The solutions 
are less trivial than those formed from Polyakov loops, since they involve kinks, but they are still 
annihilated by the volume operator. 



A somewhat different strategy was followed by Fort et al 102 |, who constructed a Hamiltonian 
lattice regularization for the calculation of certain knot invariants. They defined lattice constraint 
operators in terms of their geometric action on lattice Wilson loop states, and reproduced some 
of the formal continuum solutions to the polynomial Hamiltonian constraint of complex Ashtekar 
gravity on simple loop geometries. 



G 



2.14 The role of diffeomorphisms 

As in other discrete formulations, the spatial diffeomorphism group cannot be realized exactly on 
the lattice, and the only obvious symmetries of a cubic lattice are discrete rotations and overall 
translations. The commutator computation described above indicates that one may be able re- 
cover the diffeomorphism invariance in a suitable continuum limit. Corichi and Zapata p3| ] have 
suggested the presence of a residual diffeomorphism symmetry in the lattice theory, under which, 
for example, all non-intersecting Wilson loop lattice states would be identified. 

One can try to interpret the lattice theory as a manifestly diffeomorphism-invariant construc- 
tion, with the lattice representing an entire diffeomorphism equivalence class of lattices embedded 



in the continuum 146]. In order to make this interpretation consistent, one should modify the 
functional form of either the Hamiltonian or the measure, in such a way that the commutator of 
two lattice Hamiltonians vanishes, as a — > 0. 

2.15 The volume operator 

An important quantity in Hamiltonian lattice quantum gravity is the volume operator, the quantum 
analogue of the classical volume function J d 3 x \J det g. The continuum dreibein determinant 
detE(x) (with | detE(x)\ — det g), has a natural lattice analogue, given by 

detp(n) := ^e abc e vk p i (n,a)p : j(n,b)p k (n,c). (10) 

Apart from characterizing geometric properties of lattice quantum states, it is needed in the con- 
struction of the quantum Hamiltonian of the real connection approach. 



Loll 1 147 showed that a lattice Wilson loop state has to have intersections of valence at least 4 

in order not to be annihilated by the volume operator ^ n \l detp(n). This result is independent 
of the choice of gauge group (SU(2) or SL(2, C)). The volume operator has discrete eigenvalues, 
and part of its non-vanishing spectrum, for the simplest case of four-valcnt intersections, was first 



calculated by Loll 1 149 . These spectral calculations were confirmed by De Pietri and Rovelli [|94j , 
who derived a formula for matrix elements of the volume operator on intersections of general 
valence (as did Thiemann |188| ). 

This still leaves questions about the spectrum itself unanswered, since the eigenspaces of detp(n) 
grow rapidly, and diagonalization of the matrix representations becomes a technical problem. 
Nevertheless, one can achieve a better understanding of some general spe ctral properties of the 



lattice volume operator, using symmetry properties. It was observed in 149 that all non- vanishing 



eigenvalues of det p(n) come in pairs of opposite sign. Loll subsequently proved that this is always 



the case [ 1 5 1 1 . A related observation concerns the need for imposing an operator condition detp > 



on physical states in non-perturbative quantum gravity [151 , a condition which distinguishes its 
state space from that of a gauge theory already at a kinematical level. 

The symmetry group of the cubic three-dimensional lattice is the discrete octagonal group O, 
leaving the classical local volume function detp(n) invariant. Consequently, one can find a set 
of operators that commute among themselves and with the action of the volume operator, and 
simplify its spectral analysis by decomposing the Hilbert space into the corresponding irreducible 



representations [ 154 1 . This method is most powerful when applied to states which are themselves 
maximally symmetric under the action of 0, in which case it leads to a dramatic reduction of the 
dimension of the eigenspaces of det p{n) . 

In addition to the volume operator, one may define geometric lattice operators measuring areas 



and lengths |15CL 152 1. They are based on (non-unique) discretizations of the continuum spatial 
integrals of the square root of the determinant of the metric induced on subspaces of dimension 2 
and 1. For the case of the length operator, operator-ordering problems arise in the quantization. 

2.16 The real dynamics 

To avoid problems with the non-compactness of the gauge group and the formulation of suitable 
"quantum reality conditions" , Barbero |25|] advocated to use a real su(2)-connection formulation 
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for Lorentzian continuum gravity. This ca n be achie ved, at the price of having to deal with a more 



complicated Hamiltonian constraint. Loll |l48| , 153 translated the real connection formulation to 



the lattice and studied some of the differences that arise in comparison with the complex approach. 
Adding for generality a cosmological constant term, this leads to a lattice Hamiltonian 

tfiatt = J- 7V latt (n) [Hffi + Hpo* + VXG VdeT^)] , (11) 

n 

where schematically = (detp)- 1 / 2 ^ {Ut) P 2 , H 1 ^ = ((Tr (rUprU) - p) 2 p 6 (detp)~ 5 / 2 . This 
regularized Hamiltonian is well-defined on states with detp ^ 0, but its functional form is not 
simple. The negative powers of the determinant of the metric can be defined in terms of the 
spectral resolution of det p. The type of representation and regularization enables one to handle 
this non-polynomiality. 



2.17 Summary 

Some progress has been achieved in Hamiltonian lattice gravity, using discrete analogues of the 
Ashtekar variables. The quantization program is still open-ended, and no physically non-trivial 
solutions to the Wheeler-DeWitt equation are known. The fact that one can study and evaluate 
geometric operators provides useful characterization of quantum states. 

An analysis of the spectrum of the volume operator is crucial for handling the non-polynomial 
terms in the quantum Hamiltonian constraint. It will be necessary to find a suitable truncation or 
approximation to simplify further the spectral analysis of the Wheeler-DeWitt operator. A suitable 
quantum analogue of the continuum limit a — > has not yet been established, and in this regard 
the Hamiltonian ansatz does not go beyond the results obtained in the Lagrangian formulations 
described earlier. 

Why should one bother with a Hamiltonian quantization at all? Typical quantities one wants 
to study in a discrete path-integral approach to gravity are transition amplitudes between three- 
geometries on different spatial slices. This is not complete without a specification of the correspond- 
ing quantum states, which are in principle elements of Hilbert spaces of discrete three-geometries 
of the type described above. 



3 Quantum Regge Calculus 
3.1 Path integral for Regge calculus 



A path-integral quantization of 4d Regge calc ulus was first considered in the early eighties 181 , 104 
[Slf . This approach goes back to Regge [174], who proposed to approximate Einstein's continuum 
theory by a simplicial discretization of the metric space-time manifold and the gravitational action. 
Its local building blocks are four-simplices a . The metric tensor associated with each simplex is 
expressed as a function of the squared edge lengths I 2 of cr, which are the dynamical variables of 
this model. For introductory material on classical Regge calculus and simplicial manifolds, see 



[168, 187, 


104 


124 




197 


116 




16C, 1 


9C, 198]. 



19C, various quantum aspects are reviewed in [ 106| , 



107 



112 



110 



One may regard a Regge geometry as a special case of a continuum Riemannian manifold, a 
so-called piecewise flat manifold, with a flat metric in the interior of its 4-simplices cr, and singular 
curvature assignments to its two-simplices b (the bones or hinges). 

The Einstein action with cosmological term in the Regge approach is given by 



= J2 Vb (X — k 



A h 5. 



bObs 



V h 



(12) 



where k = g^jy, Aj, is the area of a triangular bone, Sj, = 2ir — J2az>b ^ ne deficit angle there, 

and Vb a local four-volume element. 9{o~, b) is the angle between the two 3-si mpli ces of a (it minus 
the angle between their inward normals) intersecting in b. Hartle and Sorkin [ 128 1 have generalized 



8 



( |l2| ) to the case of manifolds with boundary. The boundary contribution to the action is given by 

A#k ~ J d 3 x y/jfi) K, (13) 

frCboundary 

where ipb is the angle between the normals of the two three-simplices meeting at b. The Euclidean 
path integral on a finite simplicial complex of fixed connectivity takes the form 



Z(k,X)= [vie-^'M 



(14) 



with J VI representing the discrete analogue of the sum over all metrics. A crucial input in ( |l4| ) 
is the choice of an appropriate measure T>1. In general, a cutoff is required for both short and 
long edge lengths to make the functional integral convergent 1 104 , |8l], [l8| . One is interested in the 
behaviour of expectation values of local observables as the simplicial complex becomes large, and 
the existence of critical points and long-range correlations, in a scaling limit and as the cutoffs are 
removed. 

A first implementation of these ideas was given by Rocek and Williams [181, 183]. They 



obtained a simplicial lattice geometry by subdividing each unit cell of a hypercubic lattice into 
simplices. Their main result was to rederive the continuum free propagator (see also [101] for related 
results) in the limit of weak perturbations about flat space. This calculation can be repeated for 
Lorentzian signature [196|. 

Some non-perturbative aspects of the path integral were investigated in ]183| (see also [ 182 1 ) . 



In this work, discrete analogues of space-time diffeomorphisms arc defined as the local link length 
transformations which leave the action invariant, and go over to translations in the flat case. It 
is argued that an approximate invariance should exist in 4d. One may define analogues of local 
conformal transformations on a simplicial complex by multiplication with a positive scale factor 
at each vertex, but the global group property is incompatible with the existence of the generalized 
triangle inequalities. 



3.2 Higher-derivative terms 



Simplicial analogues of higher-derivative terms were introduced in [117, 11£]. In the continuum, 
with an appropriate choice of coupling constants, their inclusion makes the path integral less 
ill-behaved. The simplest higher-derivative term in Regge calculus is given by 



E 



A b°b 
V b 



J d 4 X V9 R^paR^, 



(15) 



with Vb denoting the local Voronoi four-volume at b. The fact that (n5| should not be identified 
with J d 4 x yfg R 2 is less surprising in the light of the classical result 81 , B2| , that Regge's expression 
Ai,Sb for the scalar curvature (for d > 2) converges to its continuum counterpart not pointwise, but 
only after integration, i.e. "in the sense of measures" (see also [103, 101], where similar convergence 
properties were studied by using an imbedding into a sufficiently large vector space IR^). 

More complicated higher-curvature terms can in principl e be constructed, using a simplicial 

pq|), but have up to now not 



analogue of the Riemann tensor (see, for example, [174, 106, 119 



been used in numerical simulations. A related proposal by Ambj0rn et al |19[| is to include terms in 
the action that depend on higher powers of the deficit angle S b , as well as terms containing powers 
of the solid angle S v at a vertex. The introduct ion of loc al vierbeins and parallel transporters is 
also necessary if one considers fermion coupling [ 104 , 170 1 . 



3.3 First simulations 

The first num erica l stu dies of the Regge action were undertaken by Berg |58|, |39| and Hamber 
and Williams [118, 11£ ]. Berg performed a Monte-Carlo simulation of the pure-curvature action 
for hypercubic 2 4 - and 3 4 -lattices with simplicial subdivision (sec [flo) for a description of the 
method) . He used the scale-invariant measure T>1 — J\ t To avoid the divergence that results 
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from rescaling the link lengths, he kept the total volume constant, by performing an overall length 
rescaling of all links after each move. This amounts to fixing a typical length scale Iq :— (vo) 1 , 
where vq is the expectation value of the 4-simplex volume. 

For k = 0, he found a negative average curvature (R), and some evidence for a canonical scaling 
behaviour of lengths, areas and volumes. For k = ±0.3, he obtained a negative (positive) average 
deficit angle ($) and a positive (negative) (R). A more detailed analysis for < k < 0.1 led Berg 
|59) to conclude that there exists a critical value k c (presumab ly a first-order transition ^(J), below 



which (R) is convergent, whereas above it diverges. (Myers [ 169 1 has conjectured that it may be 
possible to perform a similar analysis fo r M onte-Carlo data for the Lorentzian action.) 
By contrast, Hamber and Williams |118[ simulated the higher-derivative action 



S[lf] = ^ [XV b - k 5 b A b + a Ai5i/V b } (16) 

b 

on 2 4 - and 4 4 -lattices, using a time-discretized form of the Langevin evolution equation (see also 
[106 1). They employed the scale- invariant measure VI = JJ i ^F e {f), where F e enforces an ultra- 
violet cutoff e. They investigated the average curvature TZ ~ (R) and squared curvature TZ 2 (scaled 
by powers of (I 2 ) to make them dimensionless), as well as (S 2 ) and (V b )/(l 2 ). For A = k = a = 0, 
one finds a negative TZ and a large TZ 2 , indicating a rough geometry. For small a, one observes a 
sudden decrease in TZ 2 , as well as a jump from large positive to small negative values of TZ as A 
is increased. For large a, TZ is small and negative, and the geometry appears to be smooth. Like 
Berg, they advocated a fundamental-length scenario, where the dynamically determined average 
link length provides an effective UV-cutoff. 



3.4 The phase structure 

Further evidence for a transition between a region of rough and smooth geometry comes from 



Monte Carlo simulatio ns by Hamber 10£, 111| on 4 4 - and 8 4 -l attices, this time with the lattice 
VI = Yli hdh (see also [114 for a summary of results, and [|l09[ for more details on the method). 
There is a value k c at which the average curvature vanishes. For k > k c , the curvature becomes 
large and the simplices degenerate into configurations with very small volumes. He performed a 
simultaneous fit for k ci Afi and <5 in the scaling relation 

n(k) k ~ fec A n (k c - k) 5 . (17) 

This leads to a scaling exponent 6 « 0.60, with only a weak dependence on a. There are even points 
with a — that lie in the well-defined, smooth phase. Hamber also investigated the curvature 
and volume susceptibilities xn an d Xv- At a continuous phase transition, xiz should diverge, 
reflecting long-range correlations of a massless graviton excitation. The data obtained are not 
incompatible with such a scenario, but the extrapolation to the transition point k c seems somewhat 
ambiguous. On the other hand, one does not expect xv to diverge at k c , which is corroborated by 
the simulations. 



3.5 Influence of the measure 

For the pure Einstein action, (|l2] ) with A = 0, some differences between the Mi = dl/l- and the 
M2 — WZ-measure were investigated on a 4 4 -lattice by Beirl et al [3(J. A constant- volume 
constraint was used for simulations with the Mi-measure, and a cutoff I 2 < const for M^- A 
study of the ^-dependence of bulk geometric quantities agreed with previous simulations, wherever 
applicable. For small k, both measures lead essentially to identical results. For the Mi-measure 
only, one finds that (R) / (V) and (A) exhibit a jump at some fc , due to the formation of spikes 
(isolated long link lengths, with the areas staying small). From this, and the study of edge length 
distributions, one concludes that the DeWitt-type measure M2 is generally better behaved. 



3.6 Evidence for a second-order transition? 

The same W/-measure was used by Hamber in an extension of a previous simulation of the higher- 
derivative action ( |l6| ) |11S( |, for lattice sizes of up to 16 4 . At the transition point (k — k c ~ 0.244, 
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A = 1, a = 0.005), the distributions of edge lengths, volumes and curvatures are smooth and 
Gaussian-like, and the average curvature v anish es. The location of k c {a) from fits (0) to the 



average curvature coincide with those from [111], leading to 5 ~ 0.626. The data at a = do not 



seem to match this interpretation. This leads to the tentative conclusion that only for sufficiently 
large a the observed transition is of second order. It is in general "difficult to entirely exclude the 
presence of a weak first-order transition, if it has a very small latent heat". For a = 0.005, one 
finds some evidence for a decrease in the fractal dimension as k grows. 

3.7 Avoiding collapse 

Starting from the observation that there exist Regge configurations with R tot — > oo and V to t ~ * 0, 
Bcirl ct al investigated the influence of a cutoff / on the fatness of a simplex cr, defined by 

V 2 

" > / > 0. (18) 



max/ £(T / 



Such a uniform shrinking of simplices is known to be necessary in order for piecewise flat manifolds 
to approach their continuum counterparts ] jS2| . On the 4 4 -lattice, with the W/-measure and A = 0, 
they determined IZ(k) for decreasing /. For small k, the choice of / seems to have only little influ- 
ence, but towards the transition point fc c , {4> a ) simultaneously decreases. Next they investigated a 
variety of measures of the form l 2a ~ 1 dl (see also j35|). For a < 1, only a mild cr-dependence of TZ 
is observed as k — > k c . However, for a = 1.5, there are significant differences for the entire range 
of fc, and some evidence that the geometry degenerates. 

3.8 Two-point functions 

To understand the nature of the possible excitations at the phase transition, one needs to study 
correlation functions in the vicinity of k c , which is difficult numerically. Some data are available 
on the connected correlation functions of the curv atures and the volumes at fixed geodesic distance 



d, Gn(d) and Gy{d), for lattice sizes < 16 4 [115], using a scalar field propagator to determine d 



Both correlators were measured at various fc-values, leading to similar results for both a — and 
a = 0.005. The data, taken for d < 16 (~ 7 lattice spacings), can be fitted to decaying exponentials. 

Some further data (for a = 0) were reported by the Vienna group ]37|, |34|| . These authors used 
simply the lattice distance n instead of the true geodesic distance d. In |)37|, the measure was taken 
to be of the form l 2rT ~ 1 dl. They looked at Gy(n) on 3 3 x 8- and 4 3 x 16-lattices, for A = a = 1 
and A = a = 0.1, and found a fast decay for all investigated values of k, and n < 8. 

3.9 Non-hypercubic lattices 

Simulations on lattices with irregular link geometry (still with T 4 -topology) have been performed 
by Beirl et al J3l], |3^[ . They were obtained by adding a few vertices of low coordination number 
to otherwise regular lattices. One finds that the average curvature (R) increases from negative to 
positive values, even for k < k c , as a result of the formation of spikes. In [^6| , the averages (I 2 ) 
were monitored separately at the regular and the inserted vertices. 

The dependence of TZ on k is rather interesting: one observes two "critical" points, a smaller 
one k\, where the extra vertices develop spikes, and a second one ki where the remaining vertices 
follow. 1Z undergoes a small jump at k\, and a larger one at fc 2 . There is also a transition point to 
a phase with collapsed simplices at large negative fc, with a jump to large negative 1Z. Additional 
transition points at negative k were also found in simulations of the "compactified" Regge action 



S = J2b(~h Ab sm. 5 b + XVb) (this action was discussed in [ p.04| ; see also jTlJ, |140| ) and a Z2- 
version of Regge gravity ]3^ ] . Correlation functions at those points were computed in 34 1 for short 
distances, but no evidence for long-range correlations was found. 

The same authors studied the inclusion of the higher-derivative term (JlJ) in 36 . On the 
regular lattice, their findings for 7Z(k, a) confirmed those by Hamber, apart from the fact that they 
found stable expectation values even for positive TZ. Inserting irregular vertices pushes TZ to larger 
values and leads again to the appearance of an additional transition point. 
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3.10 Coupling to SU(2)-gauge fields 

Berg and collaborators ^2], |44|, |l| coupled the pure-curvature action geometrically to 

the Wilson action for SU(2)-gauge fields via dimensionless weight factors W&, 

S = S Rcggc -^Vw fc Re [Tr(l - U b )], W b = const (19) 

where Ub denotes the SU(2)-holonomy around b and /3 is proportional to the inverse square coupling 
constant, (3 = p-. One motivation was to understand whether in the well-defined pure-gravity 
region, one can choose the elementary particle masses to be <C m-pi as j3 — > oo, as one might expect 
for a realistic gravity+matter system. This seems a rather distant hope, since in the simulations 
performed so far, the ratio mpi/mhadron is of order unity. 

Initial computations were performed on a 2 x 4 3 -lattice with the scale-invariant measure, and 
at k = 0.01 and extended to larger k- values in f43| . For k < 0.04, one finds some evidence 

for a (first-order?) transition; the region of (3 where the transition occurs does not change much 
with k. Beirl et al extended this analysis by measuring the static potential V of a quark- 

antiquark pair on lattices of size 4 x 6 3 and 4 x 8 3 . With and without gravity, one finds both 
a confined and a deconfined phase; in the presence of gravity, the transition occurs at a smaller 
/3-value. More recently, Berg et al jllj have gathered further data on the location and stability of 
the well-defined phase in the (fc, /3)-plane, and extracted a string tension for various /3- values. 



3.11 Coupling to scalar matter 



The coupling of Regge gravity to scalar matter was investigated by Hamber and Williams [ 120 , 1 16| , 
for A = 1, a = 0.005, and various k- values. They considered a single, massive scalar field, with an 
action contribution 

^lE^M^yE^ (20) 

- <ij> v H i ' - i 

where Vi, Vij denote the 4- volumes associated with the vertices and edges. Most data were obtained 
on the 4 4 -lattice, with the measure y^ 1 / 30 ) Idl. The inclusion of the scalar field leads to a slight 
decrease in the average volume and edge length, as well as a very slight decrease in \1Z\. They 
again performed a fit for 1Z according to ( p7| ) and found that the critical exponent and the location 
of k c were almost unchanged. 



3.12 Recovering the Newtonian potential 

Hamber and Williams |121[ used an ansatz coming from a weak-field calculation for the amplitude 
of two gravitationally interacting particles separated by a fixed geodesic distance. They used the 
M-measure on a 16 4 -lattice, with A = 1 and a = 0. The potential V(d) is extracted from the 
connected Wilson-line to Wilson-line correlator. From data in d £ [0,5], and for sufficiently large 
k, one finds an attractive potential V{d) < 0. They considered fits to various shapes of the potential 
V, in order to extract an effective Newton's constant. - They also suggested an effective action 
for 1Z depending on k, A and a that shares some of the features observed in the simulations. 



3.13 Gauge invariance in Regge calculus? 

A view that has been expressed frequently is that away from configurations with special symmetries, 
different edge length assignments correspond to incquivalent geometries, and in this sense Regge 
calculus possesses no gauge invariance ([187, |8l], 124 1; for a dissenting opinion, see [ 184 1 ) . 

The diffcomorphism invariance can be recovered in a weak-field perturbation about flat space, 
as was shown by Rocek and Williams 1 181 , 183], and there is some evidence for the existence of 



analogous zero-modes in perturbations of regular, non-flat tesselations, at least in 2d [122]. Inspired 
by the perturbative analysis around flat space, Hamber and Williams |122| argue that a similar 
gauge invariance should persist even if one perturbs around an arbitrary non-flat background. 
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They propose as a possible definition for such gauge transformations local variations 8l 2 of the link 
lengths that leave both the local volume and the local curvature terms invariant. 

One may hope that in the non-perturbative Regge regime no gauge-fixing is necessary, since the 
contributions from zero-modes cancel out in the path-integral representation for operator averages 
[124, 122 . Mcnotti and Peirano [162, 164, 163, 161], following a strategy suggested by Jevicki and 



Ninomiya 1 1 34 ] , have argued vigorously that the functional integral should contain a non-trivial 
Faddeev-Popov determinant. Their starting point is somewhat different from that adopted in 
the path- integral simulations (see also |199| ). They treat piecewise flat spaces as special cases of 
differentiable manifolds (with singular metric), with the action of the full diffeomorphism group 
still well-defined. To arrive at a concrete representation for the Faddeev-Popov term which could 
be used in simulations seems at present out of reach. 

Recently, Hamber and Williams 1 122 , |123| have argued that the M-lattice measure is the 



essentially unique local lattice measure over squared edge lengths (this is a special case of the 
one-parameter family of local measures of the form rier[^( <T )] 1 ' 11; see a l so @] f° r a related 
derivation) . It does of course require a term with positive cosmological term in the action in order 
to suppress long edge lengths. 



3.14 Assorted topics 



Hartle [124, |125[ has suggested to compute the wave functional of the universe in a simplicial 
approximation, and to evaluate the discrete path integral semiclassically near stationary points of 
the Regge action. He investigated numerically the extrema of the action ( |l2] ) on small simplicial 
manifolds with topologies S" 4 , CP 2 and S 2 x S 2 126 1. The properti es of a Hartle-Hawking wave 
functional for a small complex with an 5 3 -boundary were studied in 127 . 



Frohlich 104] has advocated the need for a proof of reflection positivity of the Regge path 
integral, which one may expect to play a role in proving the unitarity of the theory. This can be 
formulated as a condition on the path-integral measure (including the action) under the gluing of 
two simplicial four-manifolds along a three-dimensional boundary. 

Other authors have suggested to associate gauge-theoretic instead of metric variables with the 
building blocks of a simplicial complex, for the case of the Poincare group [fnj , the Lorentz group 

and to reformulate the 



132 



[140 1, and for Ashtekar gravity with gauge group SU(2,<C) [131 
quantum theory in terms of them. 

Hamiltonian 3+1 versions of Regge calculus have been studied classically (see [190] for a re- 



view), but attempts to quantize them have not progressed very far. One meets problems with 
the definition of the constraints and the (non-)closure of their Poisson algebra. A recent proposal 
for constructing a canonical quantum theory is due to Makela [158|, who constructed a simplicial 
version of the Whcclcr-DeWitt equation, based on the use of area instead o f leng th variables (which 
however are known to be overcomplete) . In a similar vein, Khatsymovsky [141] has suggested that 
the operators measuring spatial areas ought to have a discrete spectrum. 



3.15 Summary 

Quantum Regge calculus is based on the well-explored classical discretization of the Einstein action 
due to Regge. Its weak-field limit around flat space agrees with the continuum result. Numerical 
simulations of the Euclidean path integral indicate the existence of a well-defined phase with small 
(negative) average curvature for sufficientl y small k and sufficiently large A, even in the absence 
of higher-order curvature terms. Hamber [111, 113 has found some evidence for a second-order 



phase transition in the presence of a small higher-order derivative coupling, with a vanishing average 
curvature at the transition point. These findings have not been confirmed by other groups. The 
recent controversy in the dynamical triangulations approach teaches us to treat this issue with 
some caution. 

Almost all simulations have been done on hypercubic, subdivided lattices with T 4 -topology, 
which may introduce a systematic bias in the results. There is evidence that the choice of measure 
plays a role in the appearance and suppression of singular geometries, so-called spikes. The study of 
irregular lattices suggests a direct link of the transition points with the appearance of such spikes. 
This feature is reminiscent of the appearance of singular structures in dynamical triangulations. 
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The coupling of a single scalar or SU(2)-gauge field seems to have little influence on the phase 
structure of the gravitational sector. 



4 Dynamical triangulations 
4.1 Introduction 

This quantization approach has received a lot of attention since the early nineties [l], |ll| , inspired 
by analogous studies in two-dimensional gravity, where dynamical triangulation methods have been 
a valuable tool in complementing analytical results (see, for example, |^5|, I will here exclusively 
con cent r ate on the 4d results. Other overview material is contained in Tj, 136, |195| , [|, 0, |l7], [k]. 

innm i in. 



Dynamical triangulations are a variant of quantum Regge calculus, where the dynamical vari- 
ables are not the edge lengths of a given simplicial complex, but its connectivity. A precursor is 
Weingarten's |l93| prescription for computing transition amplitudes between three-geometries, by 
summing over all interpolating four-geometries, built from equilateral 4d hypercubes living on an 
imbedding p-dimensional hypercubic lattice with lattice spacing a. Evaluating the Einstein action 
on such a configuration amounts to a simple counting of hypercubes of dimension 2 and 4, c.f. 



To avoid a potential overcounting in the usual Regge calculus, Romer and Zahringer 184 
proposed a gauge-fixing procedure for Regge geometries. They argued for an essentially unique 
association of Riemannian manifolds and equilateral triangulations that in a certain sense are best 
approximations to the continuum manifolds. The resulting "rigid Regge calculus" is essentially the 
same structure that nowadays goes by the name of "dynamical triangulations" . In this ansatz one 
studies the statistical mechanical ensemble of triangulated four-manifolds with fixed edge lengths, 
weighted by the Euclideanized Regge action, with a cosmological constant term, and optionally 
higher-derivative contributions. Each configuration represents a discrete geometry, i.e. the discrete 
version of a Riemannian four-metric modulo diffcomorphisms. At least for fixed total volume, the 
state sum converges for appropriate values of the bare coupling constants, if one restricts the 
topology (usually to that of a sphere 5 4 ). 

4.2 Path integral for dynamical triangulations 

Denoting by T the set of all triangulations of the four-sphere, the partition function for the model 
is given by 

Z(k 2 , k 4 ) = Y, cvU e ~ S[T1 ' 5 P1 = ~ K MT) + n 4 N 4 (T), (21) 
TeT * ' 

where N 2 and N4 denote the numbers of 2- and 4-simplices contained in the simplicial manifold 
T, and C(T) is the order of the automorphism group of T. One may think of ©) as a grand 
canonical ensemble, with chemical potential K4. It is related to the canonical ensemble with fixed 
volume, Z(k2, N4), by a Legendre transform 

Z( K2 ,« 4 )=^e- K ^Z( K2 ,iV 4 ) = ^ e -^ Yl e K2JV2(T) - (22) 

JV4 JV4 TeT(AT 4 ) 

The metric information is encoded in the connectivity of the simplicial decomposition, since the 
individual 4-simplices are assumed equilateral, with the edge length a set to 1. 

To understand the simple form of the action S, recall that the curvature term in Regge calculus 
(c.f.([l2"|) is represented by J^b ^5bAb, which for fixed edge length is proportional to (C4N2 — IO./V4). 
The constant c 4 = 2ir/ arccos j = 4.767 is determined from the condition that a triangulation of flat 
space should have average vanishing curvature . (Because the four-simplices a are equilateral, 

zero curvature can only be achieved upon averaging. This explains the absence of a conventional 
perturbation theory around flat space.) The cosmological term is represented by XN^Vfo) ~ XN4. 
It is sometimes convenient to re-express N 2 as a function of N , using N 2 = 2N + 2N4 — 4, valid 
for the ^-topology. The corresponding partition function is Z(kq,k 4 ) (where kq — 2k 2 ). 
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4.3 Existence of an exponential bound? 

As a consequence of identities and inequalities satisfied by the Ni (|| contains a detailed dis- 
cussion), the action (|2l]) is bounded above and below for fixed volume N4. If the number of 
configurations for fixed -/V4, is exponentially bounded as N4 — > 00, that is, Z(k2, N4) grows at most 
as Z(k2, N4) ~ e constNi^ ^gj-g j s a " cr itical line" K4 = kI(k 2 ) in the {k%, K4)-plane, where for fixed 
K2, Z{K2 1 Ki) converges for K4 > k%{k,2)- True critical behaviour may be found by approaching 
suitable points on this line from the region above the line, where Z is well defined. 

Doubts on the existence of an exponential bound were raised by Catterall et al Jf4| , who 
considered the behaviour of SI in Z(k , K4) = J2n 4 e _K4Ar4 £l( K o, N4). Their data (taken for N4 < 
32k) were consistent with a leading factorial behaviour £1 ~ (A^!)" 5 . The same scenario was favoured 
by de Bakker and Smit ]8S[] , who performed further investigations of k\. Subsequently, Ambj0rn 
and Jurkiewicz Jl2] | and Briigmann and Marinari []6Cfl added further data points at N4 = 64/c and 
N4 = 128/c respectively. Their numerical results, as well as those by Catterall et al fTq ], who 
employed an alternative method for measuring Q, favour the existence of an exponential bound, 
although they cannot claim to be conclusive. 

There have also been theoretical arguments for the existence of an exponential bound, based on 
the proofs of such bounds for the counting of minimal geodesic ball coverings of Riemannian spaces 
of bounded geometry |6f| , and the counting of discrete curvature assignments to unordered sets 
of bones || . 



4.4 Performing the state sum 

The partition function is evaluated numerically with the help of a Monte-Carlo algorithm (see j5?], 
pOj, |72fl for details). There is a set of five topology-preserving moves which change a triangulation 



locally, and which are ergodic in the grand canonical ensemble [171, |172| , 105 1 (see also ]8j for 
a discussion). No ergodic finite set of topology- and volume-preserving moves exists for generic 
four-dimensional manifolds. This prevents one from using the canonical (i.e. volume-preserving) 
ensemble. If S 4 is not algorithmically recognizable in the class of all picccwise linear (or smooth) 
4d manifolds, the numerical simulations may miss out a substantial part of the state space because 
of the absence of "computational ergodicity" |17C[ . 

An attempt was made by Ambj0rn and Jurkiewicz p3| to link the non-recognizability to the 
presence of large- N4 barriers, which should manifest themselves as an obstacle to cooling down an 
large initial random triangulation to the minimal S^-configuration. No such barriers were found 
for system sizes < 64fc, but unfortunately they were also absent for an analogous simulation (for 
N$ < 32fc) performed by de Bakker |36| for S 5 , which is not recognizable. 

Since the local moves alter the volume, one works in practice with a "quasi-canonical" ensemble, 
i.e. one uses the grand canonical ensemble Z(k 2 , K4), but adds a potential term to the action so 
that the only relevant contributions come from states in an interval [V — AN4, V + AA 4 ] around 
the target volume V. There have been several cross-checks which have found no dependence of the 
results on the width and shape of the potential term J7?| ^l) , but the lattice sizes and fluctuations 
may still be too small to detect a potential failure of ergodicity, c.f. |^3| . 

To improve the efficiency of the algorithm, Ambj0rn and Jurkiewicz fl4| used additional global 
(topology-preserving) "baby universe surgery" moves, by cutting and gluing pieces of the simplicial 
complex . In the branched polymer phase, one can estimate the entropy exponent 7, assuming a 
behaviour of the form Z(k 2 , N4) = N^ 2 ^ 3 'e K °^ K *) Ni x (1 + 0(1 IN4)), by counting baby universes 
of various sizes. At the transition point, one finds 7(^2) ~ |Io|, |lj, 0. More recently, Egawa et 
al [psi have reported a value of 7 w 0.26. 



4.5 The phase structure 

Already the first simulations by Agishtein and Migdal ^ [lj and Ambj0rn and Jurkiewicz jll| of 
dynamically triangulated 4d gravity exhibited a clear two-phase structure. After tuning to the 
infinite- volume or critical line k1(k2), one identifies two regions, K2 < and K2 > k\. The critical 
value k\ depends on the volume N4, and it was conjectured that it may move out to +00 as 
A/4 — > 00 J75|, [S9| , but it was later shown to converge to a finite value (l4|. One can characterize 
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the region with small K2 < k§ as the hot, crumpled, or condensed phase. It has small negative or 
positive curvature, large (possibly infinite) Hausdorff dimension dn and a high connectivity. By 
contrast, for K2 > one is in the cold, extended, elongated, or fluid phase. It has large positive 
curvature, with an effective tree-like branched-polymer geometry, and dn ~ 2. 

The location of the critical point on the infinite-volume line may be estimated from the peak 
in the curvature susceptibility x(«2, N A ) = ((N$) - JA 2 ) 2 )/A 4 (gj 0, [TTJ, |l9| \V2^ |(J, or the node 
susceptibility x(k ,N 4 ) = ({Nq) - (iV ) 2 )/^4 @, |§, as well as higher cumulants of N Q @. In 
Jl4) , it was suggested that one may alternatively estimate t| (-^4) by looking at the behaviour of 
the entropy exponent 7, approaching from the elongated phase. More recently, Catterall et al 
[ f77[ have used the fluctuations 
vertices as an order parameter. 



Xo = ((^0) — ( w o) j/Ni of the local volume ujq around singular 



4.6 Evidence for a second-order transition? 



It sometimes seems to be assumed that if one were to find a continuum theory at a second- 
order phase transition, it would have flat Minkowski space as its ground state (in spite of the 
S^-topology) , and gravitonic spin-2 excitations. An alternative scenario with a constant-curvature 
sphere- metric has been put forward by de Bakker and Smit ]9^, |9^[ . 

Agishtein and Migdal Q initially reported a hysteresis in the average curvature (R)(k2), indi- 
cating a first-order transition. However, subsequent authors |l91 , 192 75| , |76| found numerical 
data not incompatible with the existence of a second-order transition, and also Agishtein and 
Migdal 0] retracted their original claim as a result of a closer examination of the fixed point 
region. 

One tries to discriminate between a first- and second- (or higher-)order transition by looking 

or scaling exponents a governing the scaling behaviour 
1|, 192 1, or the peak height of susceptibilities as a function 
. Other scaling relations are discussed in 



at the Binder parameter [191, 192 
of the volume 7V4 



of suitable observables 
011 



m 73, 76 



b mm- 

However, since critical parameters are hard to measure, and it is difficult to estimate finite-size 
effects, none of the data can claim to be conclusive. 

Some doubts were thrown on the conjectured continuous nature of the phase transition by 
Bialas et al p8[ , who found an unexpected two-peak structure in the distribution of nodes near the 
fixed point. This was strengthened further by data taken at 64k by de Bakker | p7[ , with an even 
more pronounced double peak (see also J51[). Most likely previous simulations were simply too 
small to detect the true nature of the phase transition. Both Bialas et al and de Bakker observed 
that the finite size scaling exponents extracted from the node susceptibility xo grow with volume 
and may well reach the value 1 expected for a first-order transition as V — ► 00. 

The origin of this behaviour was further elucidated by Bialas et al [^9, 47 and Bialas and Burda 
[ [46| , who found a simple mean-field model that reproduces qualitatively the phase structure of 4d 
dynamically triangulated quantum gravity. With an appropriate choice of local weights, this model 
has a condensed and a fluid phase, with a first-order transition in between. A similar behaviour 
was found by Catterall et al (77), who made a related mean-field ansatz, with the local weights 
depending on the local entropies around the vertices. 



4.7 Influence of the measure 

The inclusion of a local measure term corresponding to (det g) % was studied by Briigmann and 
Marinari |5^, [57], f58|. This amounts to adding to the action a term of the form — J2 V lno(u)/5, 
where o(v) is the number of 4-simplices containing a vertex v. For /i = —5, —1, 1, 5, this seems 
to lead merely to a shift of the critical line. This setting has been revived recently by Renken 
in the context of a renormalization group analysis. 
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4.8 Higher-derivative terms 

Higher-order derivative terms were considered by Ambj0rn et al (see also |143| ), who added a 
term of the form 

to the action, where again C4 = 4.767 and 0(6) denotes the number of 4-simplices sharing a bone 
b. For h G [0,20], and with volumes up to 32fc, no major qualitative changes of the geometrical 
observables were found. The inclusion of the higher-derivative term also does not improve the 
behaviour of the average curvature (R), which continues to be positive at the critical point, whereas 
from a naive comparison with the continuum theory one would expect it to scale to zero. (This 
is also incompatible with the prediction of Antoniadis et al p0[ , should dynamical triangulations 
possess an infrared stable fixed point.) De Bakker and Smit 189] have argued that this may not be a 
reason of concern, since one expects the volume and curvature terms to mix under renormalization. 



4.9 Coupling to matter fields 

Ambj0rn et al considered the influence of both Ising spins and Gaussian scalar fields on bulk 
geometric quantities. The phase structure remained essentially unchanged, and no improvement 
in the scaling behaviour of (R) was found. Taken together with their results on higher-derivative 
gravity |^6|, one finds a universal linear dependence of the cosmological constant k^k^), with slope 
- 2.5. 

Coupling to Z 2 -spin variables s(l) located on edges I was considered by Ambj0rn et al who 
added a Wilson loop term S\v[s;T] = — (3^2 bGT o(b) Yl leb s(l) to the action. The matter sector 
behaves largely as expected when K2 is varied between the crumpled and the elongated gravity 
phase. However, in the common critical region of both sectors, where a priori one might have 
expected interesting effects, the critical behaviour seems to agree with that of the pure gravity 
system. 

More recently, Bilke et al |32| have reported a non-trivial back-reaction of matter on geometry, 
when considering coupling to several non-compact [/(l)-gauge fields. Their study was in part 
motivated by a continuum analysis of the dynamics of the conformal factor of Antoniadis et al 
[ pl| . Including three fields [/(l)-fields seems to lead to a total suppression of the branched polymer 
phase, which is replaced by a new weak-coupling phase with negative susceptibility exponent 7 
and a fractal dimension w 4. These are clearly interesting results, but should be treated with some 
caution because of the small lattice sizes involved (N4 < 16k). 



4.10 Non-spherical lattices 

Alternative topologies for the underlying simplicial complex were considered by Bilke ct al |5l], [53) . 
They showed that for topologies S 1 x S 3 and T 4 , the free energy agrees to leading order with 
that of S 4 . Neither the critical value k^(k2, N4), nor the leading contributions to the action 
density (N2)/Ni are changed. However, the simulations become more involved, since the size of 
the minimal configurations is increased. 



4.11 Singular configurations 

The singular nature of the geometry in the phase below the critical value k\ can be quantified by 
the distribution p(n) of the vertex order o(v), 

p( n ) = w(j2 s °w< n )> (24) 



first considered by Hotta et al [129]. It has a continuum part and a separate peak at high vertex 
order. One can suppress this effect by adding a term ~ J2 V (°( V ) ~~ ^) 2 to tne ac tion (see also |1| 
for a discussion of terms of a similar nature), but this leads to a simultaneous disappearance of 
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the phase transition. Hotta et al [13C] have checked that for a variety of initial configurations the 
singular structure is a generic feature of the model. 

Catterall et al observed that the pair of singular vertices form the end points of a 

singular link. They also offered a possible explanation for the formation of these singular structures: 
simplices of sufficiently low dimension can maximize their local entropy by acquiring large local 
volumes (see also p9| for a mean- field argument). Catterall et al (77) found two pseudo-critical 
points, Kq 1 ' and Kq , associated with the creation of singular vertices and links, which seem to 
merge into a single critical point reg as V — > oo. One concludes that the observed phase transition 
in the 4d dynamical triangulations model is driven by the appearance and disappearance of singular 
geometries. 

4.12 Renormalization group 

There have been attempts to apply renormalization group techniques, assumpting that the tran- 
sition is indeed continuous. Burda et al |3lJ and later Bialas et al [Q considered the "cutting of 
the last generation of minimal-neck baby universes" . This step can only be performed once, which 



severely limits the power of the procedure. Renken [178, 177 1 has applied a different blocking 
move, involving node deletion, and studied the RG flow using the volume and the vertex order as 
observables. 

4.13 Exploring geometric properties 

Since the effective geometry in both phases is rather singular, different ways of measuring length 
may lead to inequivalent definitions of "dimension" . A common (local) notion is derived from the 
volume of a geodesic ball with radius r. Usually the radius is measured in terms of the geodesic 
distance d\ (the minimal number of links), or dual geodesic distance di (the minimal number of 
links of the dual graph). Alternatively, one may consider the number n(r) of 4-simplices in spherical 
shells of thickness 1 at distance r, and define a fractal dimension by (n(r)) ~ r d F _1 jgo| [T^j 

To extract a global length scale, one may use the averages (di), (d^), or consider the average 
"radius of the universe" (r(T)) [jO] to obtain a cosmological Hausdorff dimension dcH, or the 
"average intrinsic linear extent" L = (V~ 2 J2i j di(i,j,T)} (see, for example, [[7f| [76|). Away from 
the phase transition, the fractal dimensions associated with these geometric construction are more 
or less equivalent and give dp = oo in the crumpled phase and c?f ~ 2 in the elongated phase. It 
is difficult to measure the dimension close to the transition point. 

It will not be straightforward to interpret the behaviour of observables (defined in analogy 
with the continuum theory), since in most of the phase space the geometry of the simplicial 
complex is far from approximating a metric 4-manifold. In search of a semiclassical interpretation 
for geometric observables, an alternative notion of local curvature for a simplicial manifold was 
suggested by de Bakker and Smit Q , based on a continuum expansion of the volume of a geodesic 
ball. Assuming furthermore that independent of n, n-volumes of balls with radius r behave like 
regions on S n C H n , they extracted scaling relations for various geom etric quantities for an 



intermediate range for r. This line of thought was pursued further in [185]. 

Close to the phase transition, one may investigate the behaviour of test particles (ignoring 
back-reactions on the geometry). Comparing the mass extracted from the one-particle propagator 
with the energy of the combined system obtained from the two-particle propagator |88|, ^2|, [93| , 
one does indeed find evidence for gravitational binding. 

4.14 Two-point functions 

It is possible to define two-point correlation functions on random geometries |l6|, |l4j , which are to 
be thought of as the discrete analogues of formal continuum correlators 



RicmSl 
DiffS 4 



G {r) = / Dg„ u e- S / d 4 x / d 4 x'y/det g(x)^det g(x')0(x)0(x')6(d g {x, x') - r), (25) 



for local observables 0(x), with d g denoting the geodesic distance with respect to the metric 



guv There is an ambiguity in defining the connected part of the correlator (25), as was pointed 
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out by de Bakker and Smit |9l], |9^]. Contrary to expectations, after subtraction of the square 
of the curvature expectation value, the resulting quantity (RR)(r) — (R) 2 does not scale to zero 
with large distances. They therefore proposed an alternative definition of the connected two-point 
function, by subtracting the square of a "curvature-to- nothing" correlator (R)(r). This definition 
was compared in more detail by Bialas |45| with a more conventional notion, as, for example, 
the one used in f48fl . For the case of curvature correlators, their behaviour differs significantly, 
especially at short distances. 



4.15 Summary 

In the dynamical triangulations approach, one studies the properties of a statistical ensemble 
of simplicial four-geometries a la Regge of fixed edge lengths. By summing over such discrete 



configurations according to (21), one has implicitly assumed that this leads to a uniform sampling 
of the space of smooth Riemannian manifolds. There is no obvious weak-field limit, but this is no 
obstacle in principle to the path-integral construction. Numerical simulations indicate the existence 
of a well-defined phase for sufficiently small K2 (inverse Newton's constant) and a sufficiently large 
cosmological constant. For small K2, one finds a "crumpled" phase, with small average curvature 
and a large Hausdorff dimension, and for large K2 an elongated, effectively two-dimensional polymer 
phase. At present, the consensus seems to be that the corresponding phase transition is of first 
order, with a finite average curvature at the transition point. 

Almost all simulations have been done on simplicial manifolds with 5' 4 -topology. Neither the 
inclusion of factors of (det g) p in the measure nor the addition of higher-order curvature terms to 
the action seem to have a substantial influence on the phase structure. Also matter coupling to 
spinorial and scalar fields does not seem to lead to a change of universality class, although the 
inclusion of several gauge fields may have a more drastic effect. The study of singular structures 
(vertices of high coordination number) has led to a qualitative understanding of the phase structure 
of the model. 



5 Conclusions and Outlook 

In this review, I have collected a variety of results on discrete four-dimensional models of quantum 
gravity, mainly coming from Euclidean path-integral approaches. Numerical simulations have 
yielded information on the phase structure of these models, the behaviour of two-point functions 
and a number of other properties of their partition functions. All of the path-integral models 
have some qualitative features in common. They need a (sufficiently large, positive) cosmological 
constant A to be well defined. For sufficiently small values of Newton's constant G, one finds a 
phase of collapsed geometry, with effective dimension < 4. In the gauge-theoretic model, the metric 
is degenerate, in Regge calculus, one finds spiky configurations, and in dynamical triangulations, 
the ensemble behaves like that of a branched polymer. 

In all cases, one observes a transition on the boundary of this phase, but so far no convincing 
evidence of long-range correlations has been found in its vicinity. Within the accuracy of the 
numerical simulations, this main conclusion is not altered by the inclusion of determinantal factors 
(det g) p in the measure, the inclusion of higher-order derivative terms, or the addition of matter 
fields. Why does this happen? Each of the approaches can claim that its state space represents, 
at least roughly, an approximation to the space of smooth Riemannian metrics or geometries. 
This leaves only the path-integral measure as a possible culprit. The measures used up to now 
were the simplest ones compatible with considerations of locality and gauge-invariance. It seems 
premature to blame the absence of diffeomorphism invariance (whose status in the gauge-theoretic 
formulation and the Regge calculus program remains unclear), since the explicitly diffeomorphism- 
invariant dynamical triangulations approach suffers from similar problems. Further analytical 
insights are needed to understand which modifications of the measure would make these models 
more interacting. 

There are a number of loop holes which could change the picture just presented. It is possible 
that adding enough matter of the correct type could have a non-trivial effect, or that Regge 
calculus with the inclusion of higher-order curvature terms does indeed possess a second-order 
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phase transition. Since we have very little experience with universality properties of 4d generally 
covariant theories, it is not a priori clear whether the choice of measure and the initial restrictions 
on the lattice geometry can affect the final results. 

One may of course take the attitude that something is fundamentally wrong with trying to 
construct a theory of quantum gravity via a statistical field theory approach, and that a different 
starting point is needed, an obvious candidate being a non-perturbative theory of superstrings, 
or of more general extended objects. In any case, these different approaches need not be mutally 
exclusive, and one may therefore take the results of the discrete approaches presented here as an 
indication that other attempts of constructing quantum gravity non-perturbatively may run into 
similar difficulties. 

A further unresolved problem is the "analytic continuation" of the path-integral results to 
Lorentzian signature. The Hamiltonian ansatz circumvents this problem, and some progress has 
been made in the Hamiltonian gauge-theoretic discrete approach. Although the kincmatical struc- 
ture is in place and some information on the constraint algebra has been obtained, the physical 
state space has not yet been identified. Its results are therefore not sufficiently complete to admit 
comparison with the path-integral simulations. For the simplicial formulations, only little is known 
about their canonical counterparts. One would hope that future research will throw further light 
on these issues. 
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